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Abstract 

The DFSZ model is a natural extension of the two-Higgs doublet model containing an additional 
singlet, endowed with a Peccei-Quinn symmetry, and leading to a physically acceptable axiom In 
this paper we re-examine this model in the light of some new developments. For generic cou¬ 
plings the model reproduces the minimal Standard Model showing only tiny deviations (extreme 
decoupling scenario) and all additional degrees of freedom (with the exception of the axion) are 
very heavy. Recently it has been remarked that the limit where the coupling between the singlet 
and the two doublets becomes very small is technically natural. Combining this limit with the 
requirement of exact or approximate custodial symmetry we may obtain an additional 0 + Higgs 
at the weak scale, accompanied by relatively light charged and neutral pseudoscalars. The mass 
spectrum would then resemble that of a generic two Higgs-doublet model, with naturally adjustable 
masses in spite of the large scale that the axion introduces. However the couplings are non-generic 
in this model. We use the recent constraints derived from the Higgs-WW coupling together with 
oblique corrections to constrain the model as much as possible. As an additional result, we work 
out the non-linear parametrization of the DFSZ model in the generic case where all scalars except 
the lightest Higgs and the axion have masses at or beyond the TeV scale. 
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I. INTRODUCTION 


An invisible axion M constitutes to this date a very firm candidate to provide all or 
part of the dark matter component of the cosmological budget. There are several extensions 
of the Minimal Standard Model (MSM) providing a particle with the characteristics and 
couplings of the axion [4, 5J. In our view a particularly interesting possibility is the model 
suggested by Zhitnitsky and Dine, Fischler and Srednicki (DFSZ) more than 30 years ago 
that consists in a fairly simple extension of the popular two Higgs-doublet model (2HDM). 
As a matter of fact it could probably be argued that a good motivation to consider the 2HDM 
is that it allows for the inclusion of a (nearly) invisible axion pHE] • Of course there are other 
reasons why the 2HDM should be considered as a possible extension of the MSM. Apart 
from purely aestethic reasons, it is easy to reconcile such models with existing constraints. 
They may give rise to a rich phenomenology, including possibly (but not necessarily) flavour 
changing neutral currents at some level, custodial symmetry breaking terms or even new 
sources of CP violation PHD]- 

Following the discovery of a Higgs-like particle with rrih ~ 125 GeV there have been 
a number of works considering the implications of such a finding on a generic 2HDM, to¬ 
gether with the constraints arising from the lack of detection of new scalars and from the 
electroweak precision observables m- Depending on the way that the two doublets couple 
to fermions, they are classified as type I, II or III (see e.g. P for details), with different 
implications on the flavour sector. Consideration of all the different types of 2HDM plus all 
the rich phenomenology that can be encoded in the Higgs potential leads to a wide variety 
of possibilities with different experimental implications, even after applying all the known 
phenomenological low-energy requirements. 

Requiring a Peccei-Quinn (PQ) symmetry leading to an axion does however severely 
restrict the possibilities, and this is in our view an important asset of the DFSZ model. 
This turns out to be particularly the case when one includes all the recent experimental 
information concerning the 125 GeV scalar state and its couplings. Exploring this model, 
taking into account all these constraints is the motivation for the present work. 

The structure of this paper is as follows. In section |TT| we discuss the possible global 
symmetries of the DFSZ model, namely U{1)pq (always present), and SU{2) L x SU{2) R 
(the SU(2) r subgroup may or may not be present). Symmetries are best discussed by using 
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a matrix formalism that we review and extend. Section UTT1 is devoted to the determination of 
the spectrum of the theory. We present up to four generic cases that range from the extreme 
decoupling, where the model -apart from the presence of the axiom- is indistinguishable from 
the MSM at low energies, to one where there are extra light Higgses below or around the 
TeV scale. This last case necessarily requires some couplings in the potential to be very 
small; a possibility that is nevertheless natural in a technical sense and therefore should be 
contemplated as a viable theoretical hypothesis. We discuss in detail the situation where 
custodial symmetry is exact or approximately valid in this model because the combination 
of this symmetry and naturally small couplings allows us to keep the additional scalars 
‘naturally light’ if we so wish with only one exception, meaning that the ‘contamination’ 
from the large scale present in the theory is under control. 

However, while additional scalars may exist at or just above the weak scale in this model, 


they can also be made heavy, with masses in the multi-TeV region or beyond. In section IV 
we discuss the resulting non-linear effective theory emerging in this generic situation. 

Next in section [V] we analyze the impact of the model on the (light) Higgs effective 
couplings to gauge bosons and the constraints that can be derived from the recent LHC 


data. In section [VT] we compare the potential of the DFSZ model with the most general 
potential in the 2HDM. We find out which terms of the general potential are forbidden by the 
PQ symmetry and which ones are recovered when it is spontaneously broken by the VEV 


of the scalar fields. Finally in section VII the restrictions that the electroweak precision 
parameters, particularly Ap, impose on the model are discussed. These restrictions are 
relevant only in the case where all or part of the additional spectrum of scalars is light as 
we find that they are automatically satisfied otherwise. 

We would like to emphasize that even after imposing the constraints derived from the PQ 
symmetry the model still contains enough degrees of freedom to reproduce the mass spectrum 
of a generic 2HDM, so there is no predictivity at the level of the spectrum. However, 
this nearly exhausts all freedom available, particularly if exact or approximate custodial 
symmetry is imposed. Then the scalar couplings are largely fixed and in this sense the 
model is far more predictive than a generic 2HDM —and in addition it contains the axion, 
which is its raison d’etre. 
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II. MODEL AND SYMMETRIES 


The DFSZ model contains two Higgs doublets and one complex scalar singlet, namely 



with vacuum expectation values (VEVs) («o) = ifi, (0o) = v 2 , (0) = and (aq_) = (j3 + ) = 
0. Moreover, we define the usual electroweak vacuum expectation value v = 246 GeV as 
v 2 = (y\ + v 2 )/2 and tan0 = v 2 /v\. The implementation of the PQ symmetry is only 
possible for type II models, where the Yukawa terms are 

Cy = Giq^iUR + G 2 qL(j) 2 dR + h.c., ( 2 ) 

with (f)i = ir 2 (f)*. The PQ transformation acts on the scalars as 

0i -A e* Xl 0 0 i, 02 — > e* X20 02 , 0 —> e^^cj) (3) 

and on the fermions as 

Ql Ql, II —> lLi ur —> e lXv,e uR , (1r —> e lXd6 dn, eR —> e lXed eR. (4) 

For the Yukawa terms to be PQ-invariant we need 

X u = x 1 , X d = -X 2 , X e = -X 2 . ( 5 ) 

Let us now turn to the potential involving the two doublets and the new complex singlet. 
The most general potential compatible with PQ symmetry is 

Y(0, 0!, 0 2 ) = A^(0*0 - V 2 f + A 1 (0l0 1 - Y 2 ) 2 + A 2 (0202 - Y 2 2 ) 2 

+ As(0l01 — Y ] 2 + 0202 — V ^ 2 ) 2 + A 4 (0i0l) ( 0202 ) _ (0102) (0201) 

+ ( a 0101 + ^0202)0*0 - C(0l0 2 0 2 + 0 2 010* 2 ) (6) 

The c term imposes the condition —Xi+X 2 +2 AQ = 0. If we impose that the PQ current does 
not couple to the Goldstone boson that is eaten by the Z, we also get X 1 cos 2 0 + X 2 sin 2 0 = 
0. If furthermore we choose 1 AQ = —1/2 the PQ charges of the doublets are 

Xi = — sin 2 0, X 2 = cos 2 0. (7) 

1 There is arbitrariness in this choice. This election conforms to the conventions existing in the literature. 
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Global symmetries are not very evident in the way fields are introduced above. To remedy 
this let us define the matrices [12] 


$12 — (0i 0 2 ) 


<*0 P+ 

-a_ 0o 


$21 — (02 01 ) 


00 «+ 
- 0 - «0 


= r 2 $i 2 r 2 


( 8 ) 


and 



J — $l 2 $21 — 0 2 0ll- 


(9) 


Defining also the constant matrix W = (V0 + D, 2 )I/2 + (D 2 — V 2 )t^/2^ we can write the 


potential as 


V'W, 1, J) = - y 2 ) 2 + h {Tr [(/ - W)(l + r 3 )]} 2 

+0 {Tr [(J - W)(l - r 3 )]} 2 + A 3 [Tr(/ - Wf 
+ jTr [/ 2 - (Ir 3 ) 2 ] + ill [(o + b)I + (a - b)I T3 ] 

-^Tr(J0 2 + JV* 2 )- (10) 


A SU(2)l x SU(2) r global transformation acts on our fields as 


-)• L<f>ijR\ I -)• J -)■ J. 


( 11 ) 


We now we are in a better position to discuss the global symmetries of the potential. The 
behavior of the different parameters under SU(2) R is shown in Table I. See also 0- 

Finally, let us establish the action of the PQ symmetry previously discussed in this 
parametrization. Under the PQ transformation: 

$12 -> $i 2 e ixe , 0 -> e iX * e (j) (12) 


with 


X 2 ~Xi X 2 + Xi X, A, 

A — -z-1- z -T3, Xa, — --- 


2 2 
Using the values for A 1;2 given in Eq. ([T]) 


X = 


sin 2 0 0 

0 cos 2 0 




(13) 


(14) 
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Parameter 

Custodial limit 

Ai, A 2 , A 4 

Ai = A 2 = A and A 4 = 2A 

A 3 

A 3 

ft 

^(p 

to 

v; 2 = v 2 2 = v 2 

v* 

V* 

a, b 

a = b 

c 

c 


TABLE I: In total, there are 11 parameters: 7 are custodially preserving and 4 are custodially 
breaking. See the text for our usage of the expression ‘custodial symmetry’ in the context of a 
2HDM. 

III. MASSES AND MIXINGS 


We have two doublets and a singlet, so a total of 4 + 4 + 2 = 10 spin-zero particles. 
Three particles are eaten by the W ± and Z and 7 scalars fields are left in the spectrum; 
two charged Higgs, two 0“ states and three neutral 0 + states. Our field definitions will be 


worked out in full detail in section IV Here we want only to derive the spectrum. For the 
charged Higgs mass we have at tree level 2 


m H± = 8 ( A 4 u 2 H- ±). (15) 

V S 20 J 

The quantity v$ is proportional to the axion decay constant. Its value is known to be 
very large (at least 10 7 GeV and probably substantially larger ~ 10 9 GeV if all astrophys- 
ical constraints are taken into account, see ra for several experimental and cosmological 
bounds). It does definitely make sense to organize the calculations as an expansion in v/v#. 

In the 0“ sector there are two degrees of freedom that mix with each other with a mass 
matrix which has a vanishing eigenvalue. The eigenstate with zero mass is the axion and 
A 0 is the pseudoscalar Higgs with mass 



8c 



+ V S 2 0 


(16) 


2 Here and in the following we introduce the short-hand notation s^o = sin n {mf5) and cJV = cos n (m/3). 
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Eq. (16) implies c > 0. For c = 0, the mass matrix in the O' sector has a second zero, i.e. 
in practice the A 0 held behaves as another axiom 

In the 0 + sector, there are three neutral particles that mix with each other. With hi we 
denote the corresponding 0 + mass eigenstates. The mass matrix is given in Appendix |B} 
In the limit of large v#, the mass matrix in the 0 + sector can be easily diagonalized [7] and 
presents one eigenvalue nominally of order v 2 and two of order v^. Up to 0(v 2 /v^), these 
masses are 

(ac 2 p + bs 2 p - cs 2/3 ) 2 


m 


2 hl — 32v 2 (Ai+ A 2 s^j + A 3 ) — 16u 2 


m l 2 = ~ v l + + A 2 ) - 4m 

$ 2/3 


A cj, 

: [(a - b)s 2 p + 2cc 2 pY 


.2 _ ..2 , A „.2K a - b ) S W + 2CC w} 2 , + 


m h 3 = 4Ay^ + 4u 


V — ^ C / s 29 

+ 16i> 2 


(17) 

(18) 
(19) 


A cf> 2 cj s 2 ^ A 0 

The held h\ is naturally identihed with the scalar boson of mass 125 GeV observed at the 
LHC. 


It is worth it to stress that there are several situations where the above formulae are 
non-applicable, since the nominal expansion in powers of v/v$ may fail. This may be the 
case where the coupling constants a, b , c connecting the singlet to the usual 2HDM are very 
small, of order say v/v^ or v 2 /v^. One should also pay attention to the case A^ —» 0 (we 
have termed this latter case as the ‘quasi-free singlet limit’). Leaving this last case aside, 
we have found that the above expressions for apply in the following situations: 


Case 1: The couplings a, b and c are generically of 0( 1), 


Case 2 : a, b or c are of 0(v/v^). 

Case 3: a, b or c are of (D(v 2 jvV) but c X,v 2 /v 2 . 

If c A iV 2 /v^ the 0“ state is lighter than the lightest 0 + Higgs and this case is therefore 
already phenomenologically unacceptable. The only other case that deserves a separate 
discussion is 

Case 4: Same as in case 3 but c ~ A %v 2 jv 2 ^ 

In this case, the masses in the 0 + sector read, up to 0(v 2 /v |), as 

™-l lM = 8y2 ( K T VK 2 -L ) and m 2 h3 = 4A*vJ, ( 20 ) 
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where 


A — 2 (Aic| + A 2 sjg + A 3 ) + 


cu; 


2v 2 s 


2/3 


A = 4 


Q'U 

(A 1 A 2 + A 1 A 3 + A 2 A 3 ) S' 2 ,g H- 2 (AiC^ + A 2 S^ + A 3 ) 


V 2 S2/3 


( 21 ) 


Recall that here we assume c to be of 0{v 2 /v |). Note that 


cut 


m h 1 + m /i 2 — 32u ( Aic^ + A 2 s /3 + A 3 + 


2/3 


( 22 ) 


In the ‘quasi-free singlet’ limit, when Ay —> 0 or more generically Ay -C a,b,c it is impos¬ 
sible to sustain the hierarchy v <C Uy, so again this case is phenomenologically uninteresting 
(see Appendix [C] for details). 

We note that once we set tan /3 to a fixed value, the lightest Higgs to 125 GeV and Uy to 


some large value compatible with the experimental bounds, the mass spectrum in Eq. (15), 
@ and g-g grossly depends on the parameters: c, A 4 and Ay, the latter only affecting 
the third 0 + state that is anyway very heavy and definitely out of reach of LHC experiments; 
therefore the spectrum depends on only two parameters. If case 4 is applicable, the situation 
is slightly different and an additional combination of parameters dictates the mass of the 


second (lightish) 0 + state. This can be seen in the sum rule of Eq. (22) after requiring that 
rri} l{ = 125 GeV. Actually this sum rule is also obeyed in cases 1, 2 and 3, but the r.li.s is 
dominated then by the contribution from the parameter c alone. 


A. Heavy and light states 

Here we want to discuss the spectrum of the theory according to the different scenarios 
that we have alluded to in the previous discussion. Let us remember that it is always possible 
to identify one of the Higgses as the scalar boson found at LHC, namely h\. 

Case 1. all Higgses except hi acquire a mass of order ry,. This includes the charged 
and 0“ scalars, too. We term this situation ‘extreme decoupling’. The only light states 
are hi, the gauge sector and the massless axiom This is the original DFSZ scenario 

i 

Case 2. This situation is similar to case 1 but now the typical scale of masses of 
hi, H± and A 0 is yjvv^. This range is beyond the LHC reach but it could perhaps 










be explored with an accelerator in the 100 TeV region, a possibility being currently 
debated. Again the only light particles are hi, the axion and the gauge sector. This 
possibility is natural in a technical sense as discussed in [7j as an approximate extra 
symmetry would protect the hierarchy. 


Cases 3 and 4 are phenomenologically more interesting. Here we can at last have new 
states at the weak scale. In the 0 + sector, h 3 is definitely very heavy but and 
m 2 2 are proportional to v 2 once we assume that c ~ v 2 /v 2 . Depending on the relative 
size of A i and cv^/v 2 one would have to use Eq. (17) or (20). Because in case 3 one 
assumes that cv^/v 2 is much larger than A*, h\ would still be the lightest Higgs and 
m/i 2 could easily be in the TeV region. When examining case 4 it would be convenient 


to use the sum rule (22). 


We note that in case 4 the hierarchy between the different couplings is quite marked: 
typically to be realized one needs c ~ 10” 10 Aj, where A , is a generic coupling of the 
potential. It is the smallness of this number what results in the presence of light states 
at the weak scale. For a discussion on the ‘naturalness’ of this possibility see [7j. 


B. Custodially symmetric potential 

While in the usual one doublet model, if we neglect the Yukawa couplings and set the 
I/y(l) interactions to zero, custodial symmetry is automatic, the latter is somewhat unnat¬ 
ural in 2HDM as one can write a fairly general potential. These terms are generically not 
invariant under global transformations SU(2)l x SU(2) r and therefore in the general case 
after the breaking there is no custodial symmetry to speak of. Let us consider now the case 
where a global symmetry SU(2)l x SU(2) R is nevertheless present as there are probably 
good reasons to consider this limit. We may refer somewhat improperly to this situation 
as to being ‘custodially symmetric’ although after the breaking custodial symmetry proper 
may or may not be present. If SU{2) R x SU(2) R is to be a symmetry, the parameters of 
the potential have to be set according to the custodial relations in Table [TJ Now, there are 
two possibilities to spontaneously break SU(2) L x SU(2) R and to give mass to the gauge 
bosons. 
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1. SU( 2 ) x SU( 2 ) -> 17(1) 


If the VEVs of the two Higgs helds are different (tan/3 ^ 1), the custodial symmetry 
is spontaneously broken to 17(1). In this case, one can use the minimization equations of 
Appendix [A| to eliminate V, V# and c of Eq. (10). c turns out to be of order (v/v^) 2 . In this 
case there are two extra Goldstone bosons: the charged Higgs is massless 


m H± = 0 . 


(23) 


Furthermore, the Aq is light: 


rrr Ao = 16Arr 


1 H — kS 


2 ° 2/3 


(24) 


This situation is clearly phenomenologically excluded. 


2. SU{2) l x SU(2) r -> SU(2) V 


In this case, the VEVs of the Higgs doublets are equal, so tan/3 = 1. The masses are 

m 2 H± = 8(2Xv 2 + cuj), m 2 Ao = 8 c{v 2 + uj) and m 2 h2 = m 2 H± . (25) 


These three states are parametrically heavy, but they may be light in cases 3 and 4. 
The rest of the 0 + mass matrix is 2 x 2 and has eigenvalues (up to second order in v/v, 

(a — c) 2 ' 


m hi = 16u 


A + 2 A 3 




and m 2 3 = 4 


. o A oici — c ) 2 

X^vl + Av 2 y 


*) 

(26) 


It is interesting to explore in this relatively simple case what sort of masses can be 
obtained by varying the values of the couplings in the potential (A, A 3 and c). We are 
basically interested in the possibility of obtaining a lightish spectrum (case 4 previously 
discussed) and accordingly we assume that the natural scale of c is ~ v 2 /v 2 . We have to 
require the stability of the potential discussed in Appendix [D] as well as m/ ll = 125 GeV. 
The allowed region is shown in Fig. [l] Since we are in a custodially symmetric case there 
are no further restrictions to be obtained from A p. 


C. Understanding hierarchies 

As is well known in the MSM the Higgs mass has potentially large corrections if the MSM 
is understood as an effective theory and one assumes that a larger scale must show up in the 
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Custodial Limit: Stability & M h = 125GeV 
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FIG. 1: Dark/green: allowed region in the custodial limit after requiring vacuum stability (see 
e.g. Appendix [D]) . Each point in this region corresponds to a valid set of parameters in the DFSZ 
potential. Note that c is assumed to be of order v 2 /v 2 ^ and cv^/v 2 has to be ~ Aj (case 4 discussed 
in the text). 

theory at some moment. This is the case, for instance, if neutrino masses are included via 
the see-saw mechanism, to name just a possibility. In this case to keep the 125 GeV Higgs 
light one must do some amount of fine tuning. 

In the DFSZ model such a large scale is indeed present from the outset and consequently 
one has to envisage the possibility that the mass formulae previously derived may be subject 
to large corrections due to the fact that v^ leaks in the low energy scalar spectrum. Let us 
discuss the relevance of the hierarchy problem in the different cases discussed in this section. 

In case 1 all masses in the scalar sector but the physical Higgs are heavy, of order v 
and due to the fact that the couplings A, in the potential are generic (and also the couplings 
a, b, c connecting the two doublets to the singlet) the hierarchy may affect the light Higgs 
quite severely and fine tuning of the Aj will be called for. However this fine tuning is not 
essentially different from the one commonly accepted to be necessary in the MSM to keep 
the Higgs light if a new scale is somehow present. 

In cases 3 and 4 the amount of additional fine tuning needed is none or very little. In these 
scenarios (particularly in case 4) the scalar spectrum is light, in the TeV region, and the only 
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heavy degree of freedom is contained in the modulus of the singlet. After diagonalization 
this results in a very heavy 0 + state (A 3 ), with a mass or order v^. However inspection of the 
potential reveals that this degree of freedom couples to the physical sector with an strength 
v 2 jv 2 ^ and therefore may change the tree-level masses by a contribution of order v —perfectly 
acceptable. In this sense the ‘natural’ scenario proposed in [7] does not apparently lead to 
a severe hierarchy problem in spite of the large scale involved. 

Case 2 is particularly interesting. I 11 this case the intermediate masses are of order y/vv^, 
i.e. ~ 100 TeV. There is still a very heavy mass eigenstate (/i 3 ) but again is nearly decoupled 
from the lightest Higgs as in cases 3 and 4. On the contrary, the states with masses ~ y/vv^ 
do couple to the light Higgs with strength ~ A* and thus require —thanks to the loop 
suppression factor— only a very moderate amount of fine tuning as compared to case 1 . 

It is specially relevant in the context of the hierarchy problem to consider the custodial 
case discussed in the previous section. I 11 the custodial limit the A 0 mass is protected as 
it is proportial to the extended symmetry breaking parameter c. In addition mh 2 — m#±. 
Should one wish to keep a control on radiative corrections, doing the fine tuning necessary 
to keep h\ and /12 light should suffice and in fact the contamination from the heavy fi 3 is 
limited as said above. Of course, to satisfy the present data we have to worry only about 
hi. 

IV. NON-LINEAR EFFECTIVE L AGR AN GIAN 

We have seen in the previous section that the spectrum of scalars resulting from the 
potential of the DFSZ model is generically heavy. It is somewhat difficult to have all the 
scalar masses at the weak scale, although the additional scalars can be made to have masses 
in weak scale region in case 4. The only exceptions are the three Goldstone bosons, the 
h\ Higgs and the axiom It is therefore somehow natural to use a non-linear realization 
to describe the low energy sector formed by gauge bosons (and their associated Goldstone 
bosons), the lightest Higgs 0 + state hi, and the axion. Deriving this effective action is one 
of the motivations of this work. 

To construct the effective action we have to identify the proper variables and in order 
to do so we will follow the technique described in In that paper the case of a generic 
2HDM where all scalar fields were very massive was considered. Now we have to modify the 
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method to allow for a light state (the h\) and also to include the axion degree of freedom. 
We decompose the matrix-valued $12 held introduced in Section [II] in the following form 

$12 = UM 12 . (27) 


U is a 2 x 2 matrix that contains the three Goldstone bosons associated to the breaking of 
SU{2) l (or more precisely of SU(2) L x f/(l)y to U{l) em ). We denote by G % these Goldstone 
bosons 


U = exp 


.G-t 


(28) 


Note that the matrices / and J of Eq. ([9]) entering the DFSZ potential are actually indepen¬ 
dent of U. This is immediate to see in the case of / while for J one has to use the property 
T 2 U.* = UT 2 valid for SU{2) matrices. The effective potential then does depend only on the 
degrees of freedom contained in M. l2 whereas the Goldstone bosons drop from the potential, 
since, under a global SU(2) L x SU{2) R rotation, <h 12 and U transform as 


<f> 12 ->• L^ l2 R ] U -> LUR ] =► M 12 ->• RM 12 R ] . 


(29) 


Obviously the same applies to the locally gauged subgroup. 

Let us now discuss the potential and A4 12 further. Inspection of the potential shows 
that because of the term proportional to c the phase of the singlet field <ft does not drop 
automatically from the potential and thus it cannot be immediately identified with the axion. 
In other words, the phase of the (ft held mixes with the usual 0“ scalar from the 2HDM. To 
deal with this let us find a suitable phase both in M. 12 and in (ft that drops from the effective 
potential - this will single out the massless state to be identified with the axion. 

We write M .12 — Tfi 2 U a , where U„ is a unitary matrix containing the axion. An immediate 
choice is to take the generator of U a to be the identity, which obviously can remove the phase 
of the singlet in the term in the effective potential proportional to c while leaving the other 
terms manifestly invariant. This does not exhaust all freedom however as we can include in 
the exponent of U a a term proportional to 73 . It can be seen immediately that this would 
again drop from all the terms in the effective potential, including the one proportional to c 
when taking into account that (ft is a singlet under the action of 73 that of course is nothing 
but the hypercharge generator. We will use the remaining freedom just discussed to properly 
normalize the axion and Aq fields in the kinetic terms to which we now turn. 
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The gauge invariant kinetic term will be 


/: kin = \(d„j>ydy + ^Tr [(D M $ f 12 )D^ 12 


where the covariant derivative is defined by 

12 = < 9^12 


r ) IV), ■ r'l' i+ i^rB^uTa. 


(30) 


(31) 


By defining U a — exp ( 2 ia$Xj ,/q + e-.s),j with X in Eq. (14), al! terms in the kinetic 
term are diagonal and exhibit the canonical normalization. Moreover the field a$ disappears 
from the potential. Note that the phase redefinition implied in U a exactly coincides with 


the realization of the PQ symmetry on $ 12 in Eq. (12) as is to be expected (this identifies 
uniquely the axion degree of freedom). 


Finally, the non-linear parametrization of $i 2 reads as 

$12 = UM 12 U a , 

with 

(»+ H)c„ - (S - i A 0 ) Sf V2 H +Cli 

\/ v ^+ v S 2P 

y/2H^sp (v + H)sp + (S’ + i—7 

V' 


M 1 2 — \/2 
and 


\ 


V(j> 


\ 


v 2 ,+v 2 s 2 


2/3 


-A 0 )C/3 


(32) 


(33) 


v + H = ^f[aa] + s = “ A»[a„] + %![&], H ± = C,sAt 2 , (34) 


in terms of the fields in Eq. ([Tj). The singlet field (j) is non-linearly parametrized as 


= \v c f > + p-i 


vs 2 p 


vl + V 2 s\p 


~.A 0 exp I i 


Cl(j> 


v l + v2 4p 


(35) 


With the parametrizations above the kinetic term is diagonal in terms of the fields of M 12 
and p. Moreover, the potential is independent of the axion and Goldstone bosons. All the 


fields appearing in Eqs. (33) and (35) have vanishing VEVs. 

Let us stress that H, S and p are not mass eigenstates and their relations with the hi 
mass eigenstates are defined through 


H ^ ^ A/fj/ij, S ^ ^ Rsihi, p ^ ^ Rpihi. 


(36) 


2—1 


i =1 


i— 1 


The rotation matrix R as well as the corresponding mass matrix are given in Appendix |B} 
H and S are the so called interaction eigenstates. In particular, H couples to the gauge 
fields in the same way that the usual MSM Higgs does. 
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A. Integrating out the heavy Higgs fields 


In this section we want to integrate out the heavy scalars in <3> 12 of Eq. (32) in order to 


build a low-energy effective theory at the TeV scale with an axion and a light Higgs. 

As a first step, let us imagine that all the states in $12 are heavy; upon their integration 
we will recover the Effective Chiral Lagrangian m 

2 13 

C = y Tr D^u ] D^u + a <°‘ > (37) 


i=0 


where the Oj is a set of local gauge invariant operators fT6j, and the symbol D,, represents 


the covariant derivative defined in (31). The corresponding effective couplings a* collect the 


low energy information (up to energies E ~ Anv) pertaining to the heavy states integrated 
out. In the unitarity gauge, the term D^WD^U would generate the gauge boson masses. 

If a light Higgs (h = hi) and axion are present, they have to be included explicitly as 
dynamical states mi. and the corresponding effective Lagrangian will be (gauge terms are 
omitted in the present discussion) 

~ 2 ' h h 2 


£ — ~r 1 + 2<7i— b < 72 —y + 


Tr VytfTm 


V i + V 2 S2P 


+ -dfjid^h - V ( h) 


(38) 


E 

i=0 


- 1 O t + Cr 


| | ^ 1 1 '**^ren 5 

V 


where 3 


V^U = D fl U + U(d ll U a )Ul 

formally amounting to a redefinition of the ‘right’ gauge field and 

A. 


ml 


V(h) = Y h2 ~ d:i(Xv)h 3 - d 4 -h 4 , 


jC-rcn = 77 (d ll hd l ‘h) 2 + (d.hd^h) Tr V u U ] V v U + % (d^k&'h) Tr V^VJA. 


1,4 \ p 


,2 V M 


(39) 

(40) 

(41) 


Note that the axion kinetic term is not well normalized in this expression yet. Extra contributions to the 


axion kinetic term also come from the term in the first line of Eq. (381. Only once we include these extra 


contributions, the axion kinetic term gets well normalized. See also discussion below. 


15 







Here h is the lightest 0 + mass eigenstate, with mass 125 GeV but couplings in principle differ¬ 
ent from the ones of a MSM Higgs. The terms in £ ren are required for renormalizability [H] 
at the one-loop level and play no role in the discussion. 

The couplings a* are now functions of h/v, ai(h/v), which are assumed to have a regular 
expansion and contribute to different effective vertices. Their constant parts a*(0) are related 
to the electroweak precision parameters (‘oblique corrections’). 


Let us see how the previous Lagrangian (38) can be derived. First, we integrate out from 


$12 = UMi 2 U a all heavy degrees of freedom, such as and H 0 , whereas we retain H and 
S because they contain a h\ component, namely 


$12 — !AU a Mi2, M 12 — V ^2 


(v + H)c p - Ssp 


0 


0 (v + H)sp + Sep 

where H and S stand respectively for Rhi hi and Rsi h\. 


(42) 


When the derivatives of the kinetic term of Eq. (30) act on M i 2 , we get the contribution 


dnhd^li in Eq. (38). Since the unitarity matrices, U and U a drop from the potential of 


Eq. (10) only V{h ) remains. 


To derive the first line of Eq. (38), we can use Eqs. (39) and (42) to work out from the 


kinetic term of Eq. (30) the contribution 


Tr (VJAM^V^UMa = j (1 + 2^ + .. 5 Tr VJj'V^U + C(a h h). 


(43) 


Here we used that MuM\ 2 has a piece proportional to the identity matrix and another 
proportional to r 3 that cannot contribute to the coupling with the gauge bosons since 
TrD^WD^ 1 Ut3 vanishes identically. The linear contribution in S is of this type thus de¬ 
coupling from the gauge sector and as a result only terms linear in H survive. Using that 
[U a , M 12 ] = [U a , 73] = 0, the matrix U a cancels out in all traces and the only remains of the 
axion in the low energy action is the modification D M The resulting effective action 

is invariant under global transformations U — > LUR but now R is an SU( 2 ) matrix only 
if custodial symmetry is preserved (i.e. tan/3 = 1). Otherwise the right global symmetry 
group is reduced to the U{ 1 ) subgroup. It commutes with U(1)pq. 


We then reproduce (38) with g\ = 1. However, this is true for the field H on the l.h.s. 


of Eq. (43), not h = hi and this will reflect in a reduction in the value of the gt when one 


considers the coupling to the lightest Higgs only. 
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A coupling among the S field, the axion and the neutral Goldstone or the neutral gauge 
boson survives in Eq. (43). This will be discussed in Section [VJ As for the axion kinetic 
term, it is reconstructed with the proper normalization from the first term in (30) together 


with a contribution from the ‘connection’ (d^U a )U^ in Tr T> jl U'D , ‘lA (see Eq. (52) in next 
section). There are terms involving two axions and the Higgs that are not very relevant 
phenomenologically at this point. This completes the derivation of the 0(p 2 ) terms in the 
effective Lagrangian. 

To go beyond this tree level and to determine the low energy constants aj(0) in particular 
requires a one-loop integration of the heavy degrees of freedom and matching the Green’s 
functions of the fundamental and the effective theories. 

See e.g. pans] for a classification of all possible operators appearing up to 0 (p 6 ) that 
are generated in this process. The information on physics beyond the MSM is encoded in 
the coefficients of the effective chiral Lagrangian operators. Without including the (lightest) 
Higgs field h (i.e. retaining only the constant term in the functions a,i(h/v )) and ignoring 
the axion, there are only two independent (D(p 2 ) operators 


2 2 
C 2 = l -Ti{D,UD^) +a o (0)j(Tv(T 3 U'D„U)Y 


(44) 


The first one is universal, its coefficient being fixed by the W mass. As we just saw it is 
flawlessly reproduced in the 2HDM at tree level after assuming that the additional degrees 
of freedom are heavy. Loop corrections do not modify it if v is the physical Fermi scale. The 
other one is related to the p parameter. In addition there are a few 0(p 4 ) operators with 
their corresponding coefficients 

£ 4 = ^aMgg'TiiUB^Wn - ±a s (0)g 2 Tr(UT 3 rtW^)Tr(UT 3 rfW^) + ... (45) 

In the above expression W lw and B^ u are the field strength tensors associated to the SU ( 2 ) L 
and U(l)y gauge fields, respectively. In this paper we shall only consider the self-energy, or 
oblique, corrections, which are dominant in the 2HDM model just as they are in the MSM. 

The oblique corrections are often parametrized in terms of the parameters £±, £2 and £3 


introduced in [19]. In an effective theory such as the one described by the Lagrangian (44) 


and (45) £±,£2 and £3 receive one loop (universal) contributions from the leading 0(p 2 ) term 
u 2 Tr (DJUD^W) and tree level contributions from the cij(0). Thus 


£1 — 2ao(0) + ... £2 — —g 2 a 8 (0) + ... £3 — ~g 2(l i(0) + 


(46) 
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where the dots symbolize the one-loop 0[jp 2 ) contributions. The latter are totally indepen¬ 
dent of the specific symmetry breaking sector. See e.g. [12] for more details. 

A systematic integration of the heavy degrees of freedom, including the lightest Higgs as 
external legs, would provide the dependence of the low-energy coefficient functions on h/v, 
i.e. the form of the functions a,i(h/v). However this is of no interest to us here. 

V. HIGGS AND AXION EFFECTIVE COUPLINGS 

The coupling of hi can be worked out from the one of H , which is exactly as in the MSM, 
namely 

g^HW^W* = g S i M {R H ihi + + R H3 h 3 )W li W» (47) 

where Rhi = 1 — (v/v^) 2 A\ 3 /2 and gf M = 1. With the expression of A 13 given in Ap¬ 
pendix |Bj 

9l = 9l M X ( 1 - ~^~\2 ( ac l + bs l c 2/3 - cs 2p) ) • (48) 

\ Ww / 

It is clear that in cases 1 to 3 the correction to the lightest Higgs couplings to the gauge 
bosons are very small, experimentally indistinguishable from the MSM case. In any case the 
correction is negative and g± < gf M . 

Case 4 falls in a different category. Let us remember that this case corresponds to the 
situation where c ~ \iV 2 /v 2 . Then the corresponding rotation matrix is effectively 2x2, 
with an angle 9 that is given in Appendix [B] Then 

gi — g\ M cos9. (49) 

In the custodial limit, Ai = A 2 and tan/3 = 1, this angle vanishes exactly and g 1 = g 3 AI . 
Otherwise this angle could have any value. Note however that when c A ii’ 2 /v 2 then 9 —> 0 
and the value gi ~ gf AI is recovered. This is expected as when c grows case 4 moves into 
case 3. Experimentally, from the LHC results we know [20j that g\ = [0.67,1.25] at 95% 
CL. 

Let us now discuss the Higgs-photon-photon coupling in this type of models. Let us first 
consider the contribution from gauge and scalar fields in the loop. The diagrams contributing 
to the coupling between the lightest scalar state hi and photons are exactly the same ones 
as in a generic 2HDM, via a loop of gauge bosons and one of charged Higgses. In the DFSZ 
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case the only change with respect to a generic 2HDM could be a modification in the h\WW 
(or Higgs-Goldstone bosons coupling) or in the hiH + H + tree-level couplings. The former 
has already been discussed while the triple coupling of the lightest Higgs to two charged 
Higgses gets modified in the DFSZ model to 

^fnH+H- — 8vRhi [(Ai + A 2 )s|g + ^3 + 2 A 4 ] + 16vS 2 pRsi (A 2 fig — Ais|) 

+8v (j) R p i (as 2 p + 6 c| - cs 2 p) ■ (50) 


Note that the first line involves only constants that are already present in a generic 2HDM, 
while the second one does involve the couplings a, b and c characteristic of the DFSZ model. 

The coupling of the lightest Higgs to the up and down quarks is obtained from the Yukawa 
terms in Eq. (J2]) 

C(u, d, h\) = \/2h\ [G*i ( RhiC/3 — Rsi s p ) ulUr + G 2 (RmSp + RsiCp) didn + h.c.] (51) 


The corresponding entries of the rotation matrix in the 0 + sector can be found in Ap¬ 
pendix [Bj In cases 1 , 2 and 3 the relevant entries are Rri ~ 1, Rsi ~ v2 / v p and R p i ~ v/v^, 
respectively. Therefore the second term in the first line is always negligible but the piece 
in the second one can give a sizeable contribution if c is of 0(1 ) (case 1). This case could 
therefore be excluded or confirmed from a precise determination of this coupling. In cases 
2 and 3 this effective coupling aligns itself with a generic 2HDM but with large (typically 
~ 100 TeV) or moderately large (few TeV) charged Higgs masses. 

Case 4 is slightly different again. In this case Rr i = cos6 and Rsi = sin# but R p i = 0. 
The situation is again similar to a generic 2HDM, now with masses that can be made 


relatively light, but with a mixing angle that because of the presence of the c in (B9) terms 
may differ slightly from the 2HDM. For a review of current experimental fits in 2HDM the 
interested reader can see HH. 

In this section we will also list the tree-level couplings of the axion to the light fields, 
thus completing the derivation of the effective low energy theory. The tree-level couplings 
are very few actually as the axion does not appear in the potential, and they are necessarily 
derivative in the bosonic part. From the kinetic term we get 

2R S1 


C(a 0 , h) = 




-Jiid^ap (d p G o + mzZ p ) + terms with 2 axions, 


(52) 
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From the Yukawa terms (j2j) we also get 

2 i 

C(a h q , q) = =a (j) [m u s 2 0 u^u + m d c 2 pd^d) . (53) 

^l + ^ 2 4p 

The loop-induced couplings between the axion and gauge bosons (such as the anomaly- 
induced coupling a<f,FF, of extreme importance for direct axion detection [H]) will not be 
discussed here as they are amply reported in the literature. 


VI. MATCHING THE DFSZ MODEL TO 2HDM 


The most general 2HDM potential can be read 4 e.g. from from [10| |TT 

V ( 0 i, 0 2 ) = rn 2 u 0 ^ 0 i + m 2 22 0 ^ 0 2 - m\ 2 0 j 0 2 + h.c. 

+ ( 010 l ) 2 + ( 0202 Y + A 3 ( 0101 ) ( 0 2 02 ) + A 4 | (j)\ 0 2 | 2 


Ar 


YY (0102)+ Ae (010l) (0102) + A 7 (0 2 0 2 ) (0102) + h.C 


At. 


At 


t. 


At. 


(54) 


This potential contains 4 complex and 6 real parameters (i.e. 14 real numbers). The 

most popular 2HDM is obtained by imposing a Z 2 symmetry that is softly broken; namely 
Ag = A 7 = 0 and m\ 2 ^ 0. The Z 2 approximate invariance helps remove flavour changing 
neutral current at tree-level. A special role is played by the term proportional to mi 2 . This 
term softly breaks Z 2 but is necessary to control the decoupling limit of the additional scalars 
in a 2HDM and to eventually reproduce the MSM with a light Higgs. 

In the DFSZ model discussed here v<p is very large and at low energies the additional 


singlet field 0 reduces approximately to 0 ~ v^exp^/v^). Indeed, from (35) we see that 
0 has a Aq component but it can be in practice neglected for an invisible axion since this 
component is oc v/v In addition the radial variable p can be safely integrated out. 

Thus, the low-energy effective theory defined by the DFSZ model is a particular type of 
2HDM model with the non-trivial benefit of solving the strong CP problem thanks to the 
appearance of an invisible axion 5 . Indeed DFSZ reduces at low energy to a 2HDM containing 


4 We have relabelled A —> A to avoid confusion with the potential of the DFSZ model. 

5 Recall that mass generation due to the anomalous coupling with gluons has not been considered in this 
work 
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9 parameters in practice (see below, note that v$ is used as input) instead of the 14 of the 
general 2HDM case. 

The constants A 6i7 are absent as in many Z 2 invariant 2HDM but also A 5 = 0. All these 
terms are not invariant under the Peccei-Quinn symmetry. In addition the m \2 that sofly 
breaks Z 2 and is necessary to control the decoupling to the MSM is dynamically generated 
by the PQ spontaneously symmetry breaking. There is no // = m 12 problem here concerning 
the naturalness of having non-vanishing /i. 

At the electroweak scale the DFSZ potential of eq. i§ can be matched to the 2HDM 


terms of (54) by the substitutions 


m 2 n = [—2A 1 P 1 2 + 2X 3 (V 2 + V 2 ) + anj] /4 (55) 

m 22 = [-2A ^ 2 2 + 2X 3 (V 2 + V 2 ) + bvl] /4 (56) 

ml 2 = cvl/4 (57) 

Ai = (Ai + A 3 )/ 8 , A 2 = (A 2 + A 3 )/8 (58) 

A 3 = (2A 3 + A 4 )/16, A 4 = —A 4 /16, (59) 

A 5 = 0 , A 6 = 0 , A 7 = 0 (60) 


Combinations of parameters of the DFSZ potential can be determined from the four 
masses mh x , rrih 2 , mu 0 and m#+ and the two parameters g\ (or 9) and X^h+H- that controls 
the Higgs-kPIF and (indirectly) the Higgs -77 couplings, whose expression in terms of the 
parameters of the potential has been given. As we have seen for generic couplings, all masses 
but the lightest Higgs decouple and the effective couplings take their MSM values. In the 
phenomenologically more interesting cases (cases 3 and 4) two of the remaining constants (a, 
b ) drop in practice from the low-energy predictions and the effective 2HDM corresponding 
to DFSZ depends only on 7 parameters. If in addition custodial symmetry is assumed to be 
exact or nearly exact, the relevant parameters are actually totally determined by measuring 
three masses and the two couplings (m^ turns out to be equal to m H + if custodial invariance 
holds). Therefore LHC has the potential of fully determining all the relevant parameters of 
the DFSZ model. 

Eventually the LHC and perhaps a LC will be hopefully able to assess the parameters of 
the 2HDM potential and their symmetries to check the DFSZ relations. Of course finding 
a pattern of couplings in concordance with the pattern predicted by the low energy limit 
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of DFSZ model would not yet prove the latter to be the correct microscopic theory as this 
would require measuring the axion couplings, which are not present in a 2HDM. In any case, 
it should be obvious that the effective theory of the DFSZ is significantly more restrictive 
than a general 2HDM. 

We emphasize that the above discussion refers mostly to case 4 as discussed in this work 
and it partly applies to case 3 too. Cases 1 and 2 are in practice indistinguishable from 
the MSM up to energies that are substantially larger from the ones currently accessible, 
apart from the presence of the axion itself. As we have seen, the DFSZ in this case is quite 
predictive and it does not correspond to a generic 2HDM but to one where massive scalars 
are all decoupled with the exception of the 125 GeV Higgs. 

VII. CONSTRAINTS FROM ELECTROWEAK PARAMETERS 

For the purposes of setting bounds on the masses of the new scalars in the 2HDM, £\ = A p 
is the most effective one. For this reason we will postpone the analysis of £2 and £3 to a 
future publication. 

£1 can be computed by [19] 

iW(o) n zz (o) 

£l = -14;- mT- (61) 

with the gauge boson vacuum polarization functions defined as 

n ^(q) = g^Tlvviq 2 ) + q 'V terms (V = W,Z). (62) 

We need to compute loops of the type of Fig. [2] These diagrams produce three kinds of 

X 





W / k \ V v 

AA/VWVXt AAAAAAAy 


-► \ .k + q^ / -► 

q v y q 

Y 


FIG. 2: Feynman diagram relevant for Ylyy(q). 

terms. The terms proportional to two powers of the external momentum, q^qu, do not enter 
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in Uyv(q 2 )- The terms proportional to just one power vanish upon integration. Only the 
terms proportional to k^k v survive and contribute. 

Although it is an unessential approximation, to keep formulae relatively simple we will 
compute £\ in the approximation g' = 0. The term proportional to ( g') 2 is actually the 
largest contribution in the MSM (leaving aside the breaking due to the Yukawa couplings) 
but it is only logarithmically dependent on the masses of any putative scalar state and it 
can be safely omitted for our purposes [12]. The underlying reason is that in the 2HDM 
custodial symmetry is ‘optional’ in the scalar sector and it is natural to investigate power¬ 
like contributions that would provide the strongest constraints. We obtain, in terms of the 


mass eigenstates and the rotation matrix of Eq. (|36|), 

i 9 vl 

£i = 


167t 2 u 2 
3 


^ H± 2 i 2 9 

± vi + v 2 si 


2/3 


f(m H± ,m Ao ) 


+ ( „2,\*2 S 2 /("4n m D - f( m H ± ,m 2 hi ) 

i=l \ V <P + V S 20 . 


(63) 


(64) 


where /(a, b ) = ab/(b—a ) log 6 /a and /(a, a) = a. Setting v<p —* oo and keeping Higgs masses 
fixed, we formally recover the A p expression in the 2HDM (see the Appendix in [12]), namely 

~ /K ± > m D + ^ (f( m A 0 ’ m h i ) ~ f( m H±,rn 2 h .)) 


1 


£i = 


167r 2 u 2 


2=1 


(65) 


Now, in the limit —* oo and m# ± —>• rn Ao (cases 1, 2 or 3 previously discussed) the A p 
above will go to zero as v/v$ at least and the experimental bound is fulfilled automatically. 

However, we are particularly interested in case 4 that allows for a light spectrum of new 
scalar states. We will study this in two steps. First we assume a ‘quasi-custodiah setting 
whereby we assume that custodial symmetry is broken only via the coupling A 4 b = A 4 — 2A 
being non-zero. Imposing vacuum stability (see e.g. Appendix D) and the experimental 
bound of (ei — ef M )/a = AT = 0.08(7) from the electroweak fits in [ 2 Tj one gets the 
exclusion plots shown in Fig. [3] 

It is also interesting to show (in this same ‘quasi-custodial’ limit) the range of masses 
allowed by the present constraints on AT, without any reference to the parameters in the 
potential. This is shown for two reference values of rriA 0 in Fig. [4] Note the severe constraints 
due to the requirement of vacuum stability. 

Finally let us turn to the consideration of the general case 4. We now completely give up 
custodial symmetry and hence the three masses m A0 , m H± and rrih 2 are unrelated, except 


23 













FIG. 3: Exclusion region for a ‘quasi-custodial’ 2HDM potential with a A 4 b = A 4 — 2A custodial 
breaking term. In this limit, the masses depend on the free parameters A 4 b = A 4 — 2A, A 3 and 
c = cv^/v 2 , and then the vacuum stability conditions of Appendix 
exclude regions of the free parameters above. The left, right plots are for A 3 = —0.05, A 3 = 0, 
respectively. Different color regions imply different cuts assuming that all masses (m^, m# ± and 
rrih 2 ) are greater than 100, 300 or 600 GeV (light to dark). The potential becomes unstable for 
A 3 > 0.03. 

for the eventual lack of stability of the potential. I 11 this case, the rotation R can be different 
form the identity which was the case in the ‘quasi-custodial’ scenario above. I 11 particular, 
Rs 2 = cos 9 from Appendix [B] and the angle 6 is not vanishing. However, experimentally 
cos 6 is known to be very close to one (see section [V]). If we assume that cos 9 is exactly equal 
to one, we get the exclusion/acceptance regions shown in Fig. [5j Finally, Fig. [ 6 ] depicts the 
analogous plot for cos 9 = 0.95 that is still allowed by existing constraints. We wee that the 
allowed range of masses are much more severely restricted in this case. 


0 


and AT can be used to 


VIII. CONCLUSIONS 

With the LHC experiments gathering more data, the exploration of the symmetry break¬ 
ing sector of the Standard Model will gain renewed impetus. Likewise, it is important to 
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AT= 0.08(7) at CL 95% 


AT= 0.08(7) at CL 95% 



200 400 600 800 1000 

mU GeV) 


200 400 600 800 1000 

m h ,( GeV) 


FIG. 4: Exclusion plot imposed by the constraint from AT on the second 0 + state (i.e. ‘second 
Higgs’) and the charged Higgs masses for two reference values of Ma 0 (left: 300 GeV, right: 
600 GeV) in the ‘quasi-custodial’ case discussed in Fig. |3j The concentration of points along 
approximately two axis is easy to understand after inspection of the relevant formula for AT. AT 
is vanishing at the custodial limit Mh 2 = M H ±, and also for M H ± = Ma 0 - The regions excluded 
by considerations of stability of the potential are shown. 

search for dark matter candidates as this is a degree of freedom certainly missing in the 
minimal Standard Model. An invisible axion is an interesting candidate for dark matter; 
however trying to look for direct evidence of its existence at the LHC is hopeless as it is 
extremely weakly coupled. Therefore we have to resort to less direct ways to explore this 
sector by formulating consistent models that include the axion and deriving consequences 
that could be experimentally tested. 

In this work we have explored such consequences in the DFSZ model, an extension of the 
popular 2HDM. A necessary characteristic of models with an invisible axion is the presence 
of the Peccei-Quinn symmetry. This restricts the form of the effective potential. We have 
taken into account the recent data on the Higgs mass and several effective couplings, and 
included the constraints from elcctroweak precision parameters. 

Four possible scenarios have been considered. In virtually all points of parameter space 
of the DFSZ model we do not really expect to see any relevant modifications with respect 
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AT= 0.08(7) at CL 95% 



m h2 ( GeV) 


AT= 0.08(7) at CL 95% 



m h2 ( GeV) 


FIG. 5: Exclusion plot imposed by the constraint from AT on the second 0 + state (i.e. ‘second 
Higgs’) and the charged Higgs masses for several reference values of mA 0 and tan/3 in the general 
case. The value cos 9 = 1 is assumed here. The successive horizontal bands correspond to different 
values of rriA 0 • The stability bounds have already been implemented, effectively cutting off the left 
and lower arms of the regions otherwise acceptable. 


to the minimal Standard Model predictions. The new scalars have masses of order v$ or 
^/vv^ in two of the cases discussed. The latter could perhaps be reachable with a 100 TeV 
circular collider although this is not totally guaranteed. In a third case, it would be possible 
to get scalars in the multi-TeV region, making this case testable in the future at the LHC. 
Finally, we have identified a fourth situation where a relatively light spectrum emerges. The 
last two cases correspond to a situation where the coupling between the singlet and the two 
doublets is of order v 2 i.e. very small (10 -10 or less) and in order to get a relatively 
light spectrum in addition one has to require some couplings to be commensurate (but not 
necessarily fine-tuned). 

The fact that some specific couplings are required to be very small may seem odd, but as 
it has been argued elsewhere it is technically natural, as the couplings in question do break 
some extended symmetry and are therefore protected. From this point of view these values 
are perfectly acceptable. 

The results on the scalar spectrum are derived here at tree level only and are of course 
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AT= 0.08(7) at CL 95% 



FIG. 6: Exclusion plot imposed by the constraint from AT on the second 0 + state (i.e. ‘second 
Higgs’) and the charged Higgs masses for several reference values of rriA 0 - Here we take tan/? = 1 
and allow cos# = 0.95, which is consistent with present constraints. 

subject to large radiative corrections in principle. However one should note two ingredients 
that should ameliorate the hierarchy problem. The first observation is that the mass of the 
0“ scalar is directly proportional to c; it is exactly zero if the additional symmetries discussed 
in [7J hold. It is therefore somehow protected. On the other hand custodial symmetry relates 
different masses, helping to maintain other relations. Some hierarchy problem should still 
remain but of a magnitude very similar to the one already present in the minimal Standard 
Model. 

We have imposed on the model known constraints such as the fulfilment of the bounds on 
the p-parameter. These bounds turn out to be automatically fulfilled in most of parameter 
space and become only relevant when the spectrum is light (case 4). This is particularly 
relevant as custodial symmetry is by no means automatic in the 2HDM. Somehow the 
introduction of the axion and the related Peccei-Quinn symmetry makes possible custodially 
violating consequences naturally small. We have also considered the experimental bounds on 
the Higgs-gauge bosons and Higgs-two photons couplings. Together with four scalar masses, 
these parameters determine in an almost unique way the DFSZ potential, thus showing that 
it has subtantial less room to maneuver than a generic 2HDM. 
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In conclusion, DFSZ models containing an invisible axion are natural and, in spite of 
the large scale that appears in the model to make the axion nearly invisible, there is the 
possibility that they lead to a spectrum that can be tested at the LHC. This spectrum 
is severely constrained, making it easier to prove or disprove such possibility in the near 
future. On the other hand it is perhaps more likely that the new states predicted by the 
model lie beyond the LHC range. In this situation the model hides itself by making indirect 
contributions to most observables quite small. 
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Appendix A: Minimization conditions of the potential 

The minimization conditions for the potential (|6]) are 

v 2 

Ai (2u 2 cj - V 2 ) + A 3 {2v 2 - V 2 - V 2 ) + ^ (a + ctan^) = 0, (Al) 

A 2 (2 v 2 s 2 p — V 2 ) + A 3 (2v 2 — V 2 + j- an ^^ = (A2) 

A„j ( v l ~ v l) + 2y2 ( ac i + bs l ~ cs 2/s) = 0. (A3) 

These allow us to eliminate the dimensionful parameters V\ and V 2 in favor of the 
different couplings, v and v In the case where A^ = 0 it is also possible to eliminate c. 
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Appendix B: 0 + neutral sector mass matrix 


The 3x3 mass matrix is 

( 


M HSp = 4 


8v 2 (Ai dp + \ 2 Sp + A 3 ) Av 2 (-Ai c| + A 2 sjg) s 2 p 2vv<p (acj + bs ^ - cs 2/9 ) 


\ 


4v 2 (-Ai c 2 p + A 2 s|) s 2f} ^ + 2v 2 (Ai + A 2 ) s\p -vv$ [(a - b ) s 2 p + 2 cc 2(3 ] 


* 2/3 


y 2^0 (ac 2 + bs 2 p - cs 2/3 ) -v ^ [(a - 6 ) s 2/9 + 2 cc 2/3 ] 
This is diagonalized with a rotation 


V 0 


/j/N 

w 


=i? 


//A 

h 2 

\k 3 / 


We write the rotation matrix as 


/ 

(Bl) 


(B2) 


i? = exp —A + , A T = -A, B t = -B 


and work np to second order in v/v#. We find 


A i2 — B 13 — B 23 — 0, 


so the matrix is 


R = 


i v 2 Af 3 
<P 

V 2 A.3A.3+2-Bl2 


--Alis 


u 2 A13A23— 2 j?i 2 


1 _ V 2 A 23 

vl 2 

<p 

v A 


1 - 


-Ais 
tty 1,3 

w 4 

H 2 ^13+^2: 




23 


with 


a _ 2 / 2 j t 2 \ 2 , _( a — ^) s 2/3 + 2cc 2/ g 

^13 — 7“ \ ac /3 + — CS 20) i A 23 — -2c-i-: 

A<i --A,* 


* 2/3 


(B3) 


(B4) 


(B5) 


(B 6 ) 


B \2 — - & 2 P (Aic| — ^ 2 S 2 p) + 7““X- \ 13 ( ac /3 + ^ s /3 — cs 2 ,s) [(a — b)s 2 p + 2cc 2/ g] (B7) 

C A^CZC— A^S 2/ 3 


In the case of section |I 1 1 | when the breaking of custodial symmetry is SU ( 2 ) x SU ( 2 ) — y U(l) 
the mass matrix is 

( 8v 2 [A 3 + \(sp + eg)] -2 Au 2 s 4/3 2 m^ (a + 2A ^4/s) ^ 

.2 

'2/3 


= 4 


-2An 2 s 4( g —4Au 2 cL 2A^s 4/3 


y 2vv (/) (a + 2 A^s|^ 2 A^s 4/9 


A 

'410 ' 

A«J 


(B 8 ) 


/ 
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For case 4 of section nn the rotation matrix is 

/ a a n \ 


R = 


cos 0 — sin 0 0 
sin 6 cos 6 0 
0 0 1 


tan 2 6 = 


(Al Cp - A 2 s 2 p) S 2 / 3 


(Aic| - A 2 s 2 0 ) c 2 p + A 3 - cv^/ (4n 2 s 2( 9) 


(B9) 


Appendix C: The limit A,* = 0 


The eigenvalues of the 3x3 mass matrix in the 0 + sector are 


m hi — 32?r (Aic^ + A 2 s^ + A 3 ) 

4 

2 


< = -fslpitic>+bs1) + 0(v 2 ) 


2 _ _o 2 iflC0__bSpY_ 

h3 ~ adp + bsl 


(Cl) 


Either m\ or m\ is negative. Note that the limit of a, b small can not be taken directly in 
this case. 


Appendix D: Vacuum stability conditions and mass relations 

Vacuum stability implies the following conditions on the parameters of the potential [IOj : 

Ai + A 3 > 0, X 2 + A 3 > 0, 2 A 3 + A 4 + 2a/(Ai + A 3 )(A 2 + A 3 ) > 0, 

A 3 + \J (Ai + A 3 )(A 2 + A 3 ) > 0 . (Dl) 

In the case of custodial symmetry except for A 4 B / 0, these conditions reduce to 

A + A 3 > 0, A + 2A 3 > 0, 4A + 4A 3 + A 4 b > 0 (H2) 

and assuming a, b, c very small (e.g. case 4) they impose two conditions on the masses for : 

m \ 0 + m l t ~ m l 2 > 0, m 2 H± + m 2 hl - m\ > 0. (D3) 

Appendix E: Vertices and Feynman Rules in the DFSZ model 

In the limit of g' — 0, all the diagrams involved in the calculation of £1 are of the type 
of Fig. [2] All the relevant vertices are of the type seen in Fig. [7| with all momenta assumed 
to be incoming. The relevant Feynman rules are as follows: 
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V P 

XAAAAAAx 
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FIG. 7: Two scalars and a gauge boson. 


Interaction term 

Feynman Rule for the vertex 

A V a Xd a Y 

\V a X d a Y 

A Pt- 

A (pr - PxY 


To compute II zz entering Eq. (61), we need diagrams like Fig. [2] with V 
pairs are 


Z. The X, Y 


X 

Y 

Interaction term 

Feynman Rule for the vertices 

H+ 

S 

s 

H 

G + 

C) ^ ^ tc] 

| O O I 

-\gW?H + da H. 
f /2 ? 22 W?S da A 0 

g ™^W?Sd a at 

Y V l +V S 2P 

-gW£Hd a G 0 

hgW?G + da G_ 

3 i( 2 p + q)Y 2 P + q)» 

4 ^% p {‘tp + q) lx {2p+q) v 

-g 2 0^(p + gUp + gY 

-g 2 (p + q)YP + g)v 
^(2p + q) Y 2 P + q)u 


To compute Uww entering Eq. (61), we need diagrams like Fig. [2] with V 


pairs are 


W+. The X,Y 


X 

Y 

Interaction term 

Feynman Rule for the vertices 

H + 

H + 

H + 

H 

G + 

S 

Aq 

Qj(f) 

G + 

G 0 

\gW«H_ da S 
f W+A 0 da H_- 

2 Y v i +vs 2P 

-g /r? 2 w«H_d aH 

Y V l +V S 2/3 

- gW°Hd a G _ 

yW°G 0 da G_ 

t( 2 p + q)Y 2 P + q)» 

4 v i+Xi B ( 2 P + q)»( 2 P+qY 

<t> 2 p 

g 2 ^f p (p + q)Yp + q)u 
-g 2 (p + q)YP + q)v 
t( 2 p + q) Y 2 p + q)v 
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